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Delay-independent stability in monotone systems 

Eoin Devane and loannis Lestas 


Abstract —Monotone systems comprise an important class of 
dynamical systems that are of interest both for their wide 
applicability and because of their interesting mathematical prop¬ 
erties. It is known that under the property of quasimonotonicity 
time-delayed systems become monotone, and some remarkable 
properties have been reported for such systems. These include, 
for example, the fact that for linear systems global asymptotic 
stability of the undelayed system implies global asymptotic 
stability for the delayed system under arbitrary bounded delays. 
Nevertheless, extensions to nonlinear systems have thus far relied 
primarily on the conditions of homogeneity and subhomogeneity, 
and it has been conjectured that these can be relaxed. Our aim 
in this paper is to show that this is feasible for a general class of 
nonlinear monotone systems by deriving convergence results in 
which simple properties of the undelayed system lead to delay- 
independent stability. In particular, one of our results shows 
that if the undelayed system has a convergent trajectory that 
is unbounded in all components as f — oo, then the system 
is globally asymptotically stable for arbitrary bounded time- 
varying delays. This follows from a more general result derived 
in the paper that allows to quantify delay-independent regions of 
attraction, which can be used to prove global asymptotic stability 
for various classes of systems. These also recover various known 
delay-independent stability results that are discussed within the 
paper. 

Index Terms —monotone systems; asymptotic stability; time- 
delay; nonlinear systems. 

I. Introduction 

The class of monotone systems, which comprises dynamical 
systems that preserve an order on the state space, is of signif¬ 
icant interest both due to the strong theoretical properties of 
such systems and for their suitability for modeling numerous 
physical systems. Areas in which monotonicity properties have 
frequently been exploited include population dynamics m, 
consensus protocols m, and communication systems a, 
emphasizing the breadth of applicability of monotone systems. 
The seminal papers of Hirsch, beginning with @, established a 
powerful general theory for monotone systems, demonstrating 
that the order-preservation property endows such systems with 
a rich mathematical structure. A thorough review of this theory 
can be found in 0, which details in particular the powerful 
results of generic convergence that can be obtained. 

It is well-known that delay differential equations become 
monotone under the property of quasimonotonicity 0 . In 
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the modeling of many physical problems, time-delays play an 
important role and can have a significant effect on the system 
behavior. In particular, time-delays can often destabilize a 
system and prevent convergence from occurring. It is therefore 
an important problem to ascertain how the presence of time- 
delays might affect the behavior of systems that are asymp¬ 
totically stable in the undelayed setting. Quasimonotonicity 
alone is not in general sufficient to be able to guarantee global 
asymptotic stability. Furthermore, despite the monotonicity 
property, direct comparisons between the trajectories of the 
undelayed and the delayed system are in general nontrivial to 
establish, as the presence of delays can significantly change the 
behavior of trajectories. Therefore, there have been a number 
of attempts to resolve this problem for particular types of sys¬ 
tems, through use of monotonicity combined with other system 
properties. For linear systems there exist strong results, many 
of which are based upon versions of the Perron-Frobenius 
Theorem, guaranteeing delay-independent global asymptotic 
stability and even exponential convergence rates ||2l, 0. In 
recent interesting works, these results have been extended to 
nonlinear systems satisfying properties of homogeneity a, 
in and subhomogeneity M- Furthermore, it was shown 
in ifTTll that the subhomogeneity condition can be sufficiently 
strong to imply delay-independent stability even without 
monotonicity. However, since homogeneity is closely linked to 
linearity, it is desirable to provide extensions, as was done in 
the planar case in llT2l . to a more general nonlinear monotone 
setting. 

Within this paper we investigate what conclusions can 
be drawn about the behavior of the delayed system using 
the quasimonotonicity property alone, together with simple 
properties of the undelayed system. In particular, one of our 
results is to show that whenever the system without any delay 
admits a solution that is convergent and is unbounded in all 
components in negative time, then global asymptotic stability 
is guaranteed for arbitrary bounded delays, which can be 
heterogeneous and time-varying. This follows from a more 
general result derived in the paper, which explicitly constructs 
delay-independent regions of attraction in terms of convergent 
trajectories of the corresponding undelayed system. The sig¬ 
nificance of these regions is that they allow us to generate 
proofs of global asymptotic stability for various classes of 
systems, recovering also known delay-independent stability 
results for linear/homogeneous/subhomogeneous systems and 
convergence for systems with vector fields that are negative at 
a prescribed point. 

The paper is structured as follows. The notation that will be 
used within the paper is described in Section |II] The problem 
formulation is given in Section [Till The main results are then 
stated in Section |IV] and their proofs are given in Section |V] 
Section |VT] includes a discussion on the significance of the 
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results in the paper, and conclusions are drawn in Section FVlIl 
In Appendix A we review some classical results on monotone 
systems that are used within the paper, and in Appendix B we 
derive some extensions of our results to nonpositive monotone 
systems. 

II. Notation 

The set of n-dimensional vectors with real entries is denoted 
by R", and inequalities in R" are defined as follows: x > y 
means Xi > yi for all r, x > y means Xi > yi for all i and 
X ^ y, X ^ y means Xi > yi for all i. The nonnegative 
orthant {x G R": a: > 0} is denoted by R”. 

We use C([—r, 0],R") to denote the Banach space of 
continuous functions mapping the interval [—r, 0] into R". In¬ 
equalities in C{[—r, 0], R") are treated pointwise, for example 
4> > Ip means (p{9) > ip{9) for all 9 G [—r,0]. The function 
1 represents the constant mapping 1 : [—r, 0] i-A 1„, where 
In G R" has all entries equal to 1. The product vl, where 
V G R", denotes the constant mapping vt : [—r, 0] i-A v. 
We let xt denote the function in C([—r, 0],R") given by 
Xt(9) = x{t + 9) for all 9 G [—r, 0] for a solution x{t) of the 
delayed system that will be introduced in ([2i. The subscript in 
Xk will sometimes also be used to denote the fc* component 
of the variable x{t), but this notational ambiguity will be clear 
from the context. 

Finally, the expression ^(t) will be used to denote the 
evaluation of the time-derivative of a: at f = r, ■^{x(t))\^_^. 

III. Problem formulation 

Throughout the paper we investigate the stability properties 
of a general class of monotone systems with time-varying 
delays, which we define in this section. 

We consider a function / : R” —R" and a functional 
g* : C([—r, 0],R”) ^ R" that satisfy the conditions 

a; < y and a;, = y, ^ fi{x) < /i(y), (1) 

(p <ip ^ (2) 

for all x,y € R", all <p,ip € C([—r, 0], R”), and all t > 0. 
We then formulate the delayed system 

^{t) = + g\xt), (3) 

which will be analyzed throughout the paper. We also consider 
the corresponding undelayed system 

n‘T 

-{t) = f{x{t))+g\x{t)t). ( 4 ) 

We assume that the time-dependence t i-a g^{4>) is continu¬ 
ous for all hxed p G C([—r, 0], R") and that the mappings 
X I— f{x), (p I— g^ip) are both locally Lipschitz for all 
hxed t > 0. The systems (|2l and © are then guaranteed 
to have unique solutions throughout R+ x R" which exhibit 
continuous dependence on the initial data lfT3l Theorems 2.2.1, 
2.2.2, and 2.2.3]. 

It is well known that conditions O and © render sys- 
teirQ © monotone, i.e. if two initial conditions Xta^Xta G 

^Note that system E) is also monotone as this is a special case of 0. 


C([—r, 0],R") satisfy Xt^ < Xtg, then the corresponding 
solutions of © satisfy x{t) < x{t) for all f > fo ©■ This 
follows from the fact that the vector held in © satishes the 
property of quasimonotonicity (see more detailed discussion 
of this in Appendix A). 

Furthermore, it is assumed that given any v G R", the map 
1 1 —> y*(ul) is constant, independent of time t. This is always 
satished if the time-variation in y* is due to time-varying 
delays. Therefore, system ©, whose stability properties will 
be investigated, is a monotone system that can have arbitrary 
bounded time-varying delays. 

The equilibria of the systems © and @ are known to 
be equivalent 0. The problem we are interested in is to 
relate the stability properties of an equilibrium of the delayed 
system ©, which is inhnite-dimensional, to those of this 
equilibrium for the hnite-dimensional undelayed system ©. 
For simplicity, we will suppose that the equilibrium under 
consideration is located at the origin, whence © and @ 
become positive systems (i.e. the nonnegative orthant R" is 
positively invariant with respect to both systems) due to the 
monotonicity properties © and ©. As is common in the 
analysis of positive systems, we use standard definitions for 
stability (e.g. Definition 5.1.1 in lfT3l l but with any region for 
Xt used in the definition restricted to the cone C{[—r, 0], R"). 
The origin is defined to be a globally asymptotically stable 
equilibrium of © if it is stable and limt_i.oo x{t) = 0 for any 
solution xit) of © with initial condition xt^ in C{[—r, 0], R"). 
In Appendix B we show how the stability results in Section HVl 
can then be extended to establish convergence results also for 
monotone systems that are not positive. 

IV. Convergence results 

In this section we will state the main results of the paper. 
The proofs of these results are given in Section |V] and further 
discussion on their significance is included in Section |VT] 

One of our main results is given in the following theorem, 
which states that the origin in © is globally asymptotically 
stable if there exists any single convergent trajectory of the 
undelayed system @ that is unbounded in all components in 
negative time. 

Theorem 1: Suppose that the undelayed system @ ad¬ 
mits a solution y{t) satisfying limt_>._oo = oo and 

limt_>.oo = 0 for all components i. Then the origin is 
globally asymptotically stable for the delayed system ©. 

Remark 1: It should be noted that system © is allowed to 
have arbitrary bounded delays, which can be heterogeneous 
and time-varying. 

Remark 2: Theorem[T]can be seen as addressing the problem 
of when global asymptotic stability of an undelayed mono¬ 
tone system also implies delay-independent global asymptotic 
stability. As discussed in the introduction, existing results 
in the literature have relied primarily on the properties of 
homogeneity and subhomogeneity, and it has been conjectured 

^By positive invariance of a subset X C R'* with respect to 0, we mean 
that Xfp g C([—r, 0],JV) implies x{t) G X for all t > Iq. Likewise, when 
we say that Y C R" is a region of attraction of an equilibrium x* for 0, 
we mean that xt^ € C([—r, 0], Y) implies limt_>.oo x{t) = x*. 
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in m that these assumptions can potentially be relaxed. 
Theorem [T] shows that this is possible for systems having a 
backward-time divergent trajectory. 

Remark 3: The conditions on y{t) given in Theorem [T] 
trivially imply global asymptotic stability for an undelayed 
monotone system. Nevertheless, proving that stability is not 
compromised under arbitrary bounded delays is a more in¬ 
volved problem. This is due to the fact that direct comparisons 
between the trajectories of the undelayed and the correspond¬ 
ing delayed system are in general nontrivial to establish. A 
main part of the proof of the theorem is therefore associated 
with the derivation of results that allow us to quantify such 
connections. 

Theorem [T] follows from Lemma [T] stated below, which 
is a more general technical result that allows the explicit 
deduction of delay-independent regions of attraction for Q 
from convergent trajectories of the undelayed system Q. 

Lemma 1: Suppose that the undelayed system (01 admits a 
solution y{t) satisfying yi{0) > 0 and limj^oo 2/i(i) = 0 for 
all components i. Let > 0 be a point defined in terms of y{t) 
by means of equations (|5]l-(l8]). Then any solution of the de¬ 
layed system ([3]) with initial condition Xtg € C([—r, 0], [0,^]) 
satisfies limt_>oo x{t) = 0. 

Remark 4: The significance of the point C, will become 
clearer in the proof of Lemma [T] and in the discussion in 
Section |Vl] In particular, a main significance of Lemma [T] 
is that for many classes of systems the point C, can be chosen 
to be arbitrarily large, thereby guaranteeing delay-independent 
global asymptotic stability. This leads to the proof of Theo¬ 
rem [T] and it also recovers known delay-independent stability 
results for linear/homogeneous/subhomogeneous systems. 

Lemma[T]has the following immediate corollary for systems 
with continuously differentiable vector fields. In particular, 
it is guaranteed in this case that ^ 0, which ensures 
delay-independent asymptotic stability of the system (O. Note 
that this condition on the vector field is not necessary to 
deduce asymptotic stability and is also not an assumption in 
Theorem [T] 

Corollary 1: Suppose that / and are both continu¬ 
ously differentiable for all fixed t and that the undelayed 
system 01l admits a solution y(t) satisfying yi(0) > 0 and 
yi{t) = 0 for all components i. Then the point ^ 
in Lemma [1] can be chosen to be strictly positive in all 
components, and the origin is asymptotically stable for the 
delayed system ([3). 

Corollary [T] implies also the following simple corollary. 

Corollary 2: Suppose that / and are both continuously 
differentiable for all fixed t and that the origin is asymptoti¬ 
cally stable for the undelayed system (0|i. Then the origin is 
also asymptotically stable for the delayed system (O. 

V. Proofs of the results 

We now present the proofs of the main results stated in 
Section EYl We begin with the derivation of Lemma [T] and 
then use this to prove Corollaries [T] and |2] and Theorem [T] 

Within the proof of Lemma [T] we will make frequent use 
of two important results from ||6l, namely Theorem 5.1.1 and 


Corollary 5.2.2. These results are quoted in the appendix as 
Lemmas [3and|4l 

Proof of Lemma [7} We first define, in terms of the 
known trajectory y{t) of (01, various quantities that will be 
used throughout the proof of the lemma, and which are also 
needed in order to define the point Q. 

• Firstly, for all x < j/(0), we define, as in lfT4l . the 
functions 

Ti{xi) := sup{r: Xi < y^{s) \fs G [0,t]}, 

Vi{x^) := (5) 

T{x) := mhiTi{xi), and V{x) := uiaxVi{xi) 

i i 

in terms of the solution y{t) of (0]). Then V{x) = 
and, according to the given asymptotic properties of y{t), 
V : [0,2/(0)] —>■ IR,+ is well-defined. Moreover, due to the 
continuity of trajectories of (|4|l, each Ti, and hence Vi, 
is left-continuous. 

« We also introduce the mapping hi : R+ —R+, 

:= infjr > u: 3a G [0,2/i(0)] 

(o) 

such that Ti{a) = r}, 

which is well-defined by the asymptotic proper¬ 
ties of y{f). 

• Let fP denote a time for which 

2/(fJ’) <2/(0). (7) 

Using the map h in (|6]), we define the point G R" by 

Q ■■= y^{h^{n) ( 8 ) 

for all i. Not410 that C > 0. 

The significance of these notions will become clear throughout 
the proof. We now give a brief sketch of the proof to facilitate 
its readability. The function V{x) will be used to construct 
a compact positively invariant set S for the delayed system. 
We will then show that is a maximal point in this set. 
This property will be used to show that the vector field is 
nonpositive at C,, for an appropriately constructed bounding 
monotone system. Hence the desired convergence result can 
then be deduced. 

The proof is now given in several steps: 

Step 1: We first show that the level set 

5:={a;G [0,2/(0)]: U(a:) <e-*"} (9) 

is positively invariant with respect to the delayed system @- 

To establish this, we use an argument similar to Lyapunov- 
Razumikhin analysis to demonstrate that any trajectory ini¬ 
tially within S must remain within 5 for all future time. We 
begin by proving the following lemma. 

^If the vector field is also continuously differentiable, then it follows from 
Corollary ^ that can be chosen such that 0. It can easily be seen from 
the definition of C, that the weaker condition y{t) 0 Vt C [0,f^] is also 
sufficient for the existence of a ^ that is strictly positive. 

"^It should be noted that for a given trajectory y{t) there could exist multiple 
points (and times tP) that satisfy (7), with each such point leading to 
a different estimate of the region of attraction. It will be seen in the proof 
of Theorem ^ in Corollary |4| and in the example in Section IVI-CI that for 
various classes of systems the point can be appropriately chosen so as to 
also deduce global asymptotic stability. 
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Lemma 2: Suppose that y(a;(s + 0)) < y(a:(s)) for all 
9 S [—r, 0], for a solution x{t) of the delayed system Q and 
a time s for which the segment Xg € C{\—r^ 0], [0, j/(0)]). Let k 
be any component of x{s) such that F(x(s)) = Vk{xk{s)). 
Then 


Proof of Lemma |2} The conditions stated in Lemma |2] 
imply that T{x{s + 6)) > T{x{s)) for all 9 G [—r, 0]. This 
means that whenever x{s) < y{a) for all a € [ 0 , s], we must 
have x{s + 9) < y{s). In particular, for s = T{x{s)) we 
therefore get x{s + 9) < y{T{x{s))) for all 9 G [—r, 0]. Now, 
by definition, y(a:(s)) = Vk{xk{s)) implies yk{T{x{s))) = 
Xk{s), whence ([T]i and (|2l) respectively give 

fk{xis)) < fk{y{T{x{s)))), (10) 

gkixs) < gi{yiT{x{s)))t). (11) 

By the particular choice of component k and the definition (|5]l, 
we have ^^(T(x(s))) < 0, whence the properties (fTOl i 
and (fTTl i and the time-invariance of g* for constant arguments 
give 

^(s) < fk{yiTix{s)))) + gl{y{T{x{s)))l) 

= fk {y(T{x{s)))) + Pfc (y(r(x(s)))l) 

= ^(r(x(s))) < 0 . 


Let us now consider the set defined in (|9]l. The asymptotic 
properties of y{t) and the left-continuity of V ensure that S 
is compact. By this compactness, the only way in which the 
desired positive invariance could be violated is if there exists 
a trajectory x{t) of OJ and a time r such that x{a) G S 
for all (T < T and for which given any e > 0 there exists 
V G {t,t + e) with x{v) ^ S. Since V{x{t)) < e~*’’ implies 
T{x{t)) > tP and we know that yit^) <gC y(0), it follows 
that x{t) <C y(0). We can thus always choose e sufficiently 
small that x{a) G [0,j/(0)], so that V{x{a)) is defined, for all 
(T G [r, u]. Therefore it must be the case that V(x(cr)) < 
for all (T < T and V(x(v)) > e From the continuous 
differentiability of the trajectory x(t) and the left-continuity 
and strictly increasing nature of each Vi, it follows that there 
must exist an interval ( 71 , 72 ) C [t,v\ such that V{x{'-fi)) > 
max{e“*’’,sup^g[,-..y^] F(x(ct))} and on which the mapping 
a I—>■ V{x{a)) is strictly increasing. Moreover, we know from 
the continuous differentiability of the trajectory x = y{t) of (|4|i 
that each Vt is continuous everywhere except on a (possibly 
infinite) discrete set of isolated points. This means that there 
must exist some 73 G ( 71 , 72 ) such that l/(x(cr)) = Vk{xk{<j)) 
for all cr G ( 71 , 73 ) for some fixed choice of k. It thus follows 
that the mapping tr 1 —>■ Xk{cr) is strictly increasing on ( 71 , 73 ), 
whence the Mean Value Theorem ensures the existence of 
some /3 G ( 71 , 73 ) at which Xk(/3) > 0. But, by construction, 
we know that V{x{f})) > V{x{f} -\- 9)) for all 9 G [—f, 0] 
and V(x(/3)) = 14(x/c(/3)), so this positive time-derivative is 
in contradiction to Lemma |2] Therefore, S must be positively 
invariant with respect to Q. 


Step 2: We now show that is a maximal point of the set 
S. In particular, we show that S C [0, and ( G S. 

To do this, consider any x G S. According to the defini¬ 
tion (| 6 ]l, it must follow from Ti{xi) > that Ti{xi) > hi{t'P), 
whence Xi < yi{hi{tP)) = (i. Since this holds for all i, we 
immediately conclude that S C [0, C]- 

Additionally, since each T) is a nonincreasing, left- 
continuous function on a compact set, the image Ti([ 0 , t/i( 0 )]) 
contains all of its right-limits. Therefore, by ©, for any i 
there must exist some a G [0, yi(0)] such that Ti{a) = hi{tP). 
It thus holds by (|5]l that a = yi{hi{tP)) = Q. This shows that 

lies within the domain of definition of T, V and moreover, 
because hi{tP) > fP, that T(Q > fP. Therefore, G S. 

Step 3: We now use the results of Steps 1 and 2 to guarantee 
that the set [0, C] lies within the region of attraction of the 
origin of 0 . 

From Steps 1 and 2, we conclude that the sublevel set S 
satisfies positive invariance with respect to Q, is contained 
within the interval [0, (^], and contains the point (. Applying 
positive invariance to the trajectory of Q with initial condi¬ 
tion (jl at time 0 immediately implies that the vector field at 
time 0 evaluated at ^ must satisfy 

/(C)+5°(C1) <0. ( 12 ) 

To make use of (fT 2 l i. we define the functional 

5 :C([-r, 0 ],R")^lR”, 

gW'=g°{ sup (I>i9)l), (13) 

SGl-r-.O] 

where the supremum is taken componentwise, i.e. 

supgg[_^ Q] ())(0) is a vector in R" with components 
suP 6 (ie[_r 0 ] ^^i(^j)- *0 is time-invariant, locally Lipschitz, 
and satisfies the order-preserving property (|2]l. Thus we can 
form the time-independent monotone system 

^{t) = f{x{t))+g{xt), (14) 

which has an equilibrium at 0 and is hence positive. More¬ 
over, (Ell shows that /(C)+ 5 (Cii) — 0 ^ whence an application 
of Lemma |4] guarantees that the solution z{t) of (fT4l) with 
initial condition zt^ = Cl at any initial time to must satisf}@ 
limt^oo z{t) = 0. 

We then observe that, by © and the time-invariance of 
5 * for constant arguments, it must always hold that g*{(j>) < 
g*(supeg[_^^o] f{9)t) = 5 °(supeg[_^_o] ^{9)1) = gif). Thus, 
the vector fields of the monotone systems © and Gli are 
related by 

fifio))+g\f)<fim)+m (15) 

for all f > 0 and all f G C([—r, 0], R”). Inequality ( fTSl ) 
means that we may invoke Lemma [3 to deduce that if xit) 

^To be precise, since positivity means that the trajectory is kno\vn to be 
bounded belo\v by 0 in all components, Lemma[4]gives limt —>.00 z(t) = x* 
for some x* satisfying 0 < x* < z{to) = (^ <C y(0). The continuity of the 
maps / and g then means that any such point must satisfy f(x*)-\-g{x*l.) = 
0 and so f(x*) -h g^(x*l.) = 0, meaning that it is an equilibrium of the 
undelayed system s. But then, because y(0) ^ 0 and limt^oo y(i) = 0, 
the monotonicity of the system 0 prohibits this equilibrium from lying in 
0 < X* < y(0), whence we conclude that x* = 0 is the only possibility. 
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denotes the solution of @ through any xtj, G C{[—r, 0], [0, C]), 
then it is guaranteed that x{t) < z{t) for all t > to. 
But since lirut^oo z{t) = 0 and x{t) must always remain 
nonnegative, this immediately implies that limt_).oo x{t) = 0. 
This proves that all trajectories with initial conditions within 
the set C{[—r, 0], [0, C]) converge to the origin. ■ 

We now prove Corollary[Tl which states that when the vector 
field is continuously differentiable then the region of attraction 
in Lemma [T] has a nonempty interior and asymptotic stability 
can be deduced. 

Proof of Corollary |7} Since / and p* are continuously 
differentiable, 0 Remark 3.1.2] guarantees that y{t) ^ 0 
for all t > 0. Let be a time for which y(t^) ^ 2/(0). 
Since y{t) » 0 for all t G [0,f^] and limt_>.oo 2/(f) = 0, 
there must exist t > such that Ti{a) = r for the choice 
a = ^ inftg[o,tJ>] yi{t) > 0. This holds for all i, so we see that 
the point defined in ® satisfies ^ 2> 0. 

In addition, stability follows from convergence, using the 
monotonicity of the system Q and the continuous dependence 
of its solutions on their initial conditions. We thus see that 
the origin is both stable and attractive on the set [0, C], 
where ^ 0. We therefore conclude that the origin is an 
asymptotically stable equilibrium for the delayed system Q 
with arbitrary bounded delays, with a region of attraction 
containing the nonempty interval [0, ()]. ■ 

Corollary |2] now follows directly from Corollary [T] 

Proof of Corollary \2}; Asymptotically stability of (|4]i 
means that there exists e > 0 such that for all y(0) G (0, e), 
the solution y{t) of (|4|l satisfies limt_>oo y(t) = 0. Corollary[T] 
thus applies. ■ 

We finally show how the method employed in proving 
Lemma [U enables us to prove the global stability result stated 
in Theorem [T| Note that continuous differentiability is not 
required here. 

Proof of Theorem Q} Consider the modified definitions 
Ti{xi) := sup{r: Xi < yi{s) Vs G (-oo,t]}, 

V,{x,) := 

T{x) := minTi(xi), and V{x) = maxl^(xi), 

i i 

in terms of the solution y{t) of (|4|i. Then V{x) = 
as before and V is now a well-defined, positive defi¬ 
nite, left-continuous function on the whole of IR", accord¬ 
ing to the given asymptotic properties of y{t) and the conti¬ 
nuity of the system © 03 . The structure of the proof is then 
analogous to that of Lemma [T] 

Step 1: As a consequence of the fact that V is now de¬ 
fined on the whole positively invariant orthant R", analogous 
arguments to those in Step 1 of the proof of Lemma [T] show 
that all sublevel sets Sc ■= {x G R" : V(x) < c} are now 
positively invariant with respect to ©. 

Step 2: Consider an arbitrary element x G Sc- Following 
the proof of Lemma [T] we introduce the mapping hi{u) := 
inf{r > u: 3a G R+ such that Ti(a) = r}. Then, for any 
X G Sc we have V{x) < c, implying that Ti{xi) > —logc 
for all i. Thus, it must hold that Ti{xi) > hi{— logc), and so 
Xi < yi{hi{ —logc)). Consequently, if we define Cf G R" by 

Q ■=yi{ht{-logc)) (16) 


for all L then Sc C [ 0 ,^°]. 

In particular, the sublevel sets are thus bounded, meaning 
that we can again use the left-continuity and nonincreasing 
nature of Ti to deduce the existence of a G R+ such that 
hi{—\ogc) = Ti{a). Therefore, a = yi{hi{—logo)), whence 
we have all Ti{Q) = hi{— logc) > — logc and so V{Cf) < c. 
This implies that G Sc- 

Step 3: As in the proof of Lemma [T] we therefore have 
an invariant set Sc containing and contained within [ 0 , C°], 
whence we deduce 

/(C)+5°(ri) <0. (17) 

Furthermore, we observe that since limi_i._oo yi(f) = cjo 
for all i, the function V is radially unbounded. Therefore, 
the sublevel sets satisfy Sc — >• R+ as c —> oo. In particular, 
the points cxd as c —> cxd, meaning that there exist 

arbitrarily large for which dnii holds. Thus, given any 
Xto G C([—r, 0], R"), there exists such that xtg < 
and (fTTl l holds. Then repeating the construction of the bound¬ 
ing system (fl4ll and applying Lemmas |3] and @1 as in Step 
3 of the proof of Lemma [T] guarantees that the solution xif) 
of © with initial condition xt^ must satisfy limi_>oo x{f) = 0 . 
Stability now follows by the argument used in the proof of 
Corollary [T] Therefore, the origin is globally asymptotically 
stable for the delayed system ©. ■ 

Remark 5: It should be noted that even though the condi¬ 
tion limt_>_oo yi(f) = oo for all components i immediately 
implies the existence of arbitrarily large points at which each 
component of the vector field is individually nonpositive, it 
is nontrivial to show that this happens simultaneously for all 
components. The latter is needed in order to deduce delay- 
independent stability using Lemma |4] and the monotonicity 
arguments in Lemma |3] It was shown above that the quantity 
(f constructed within the proof of Theorem[T] explicitly defines 
such points, thus allowing global asymptotic stability to be 
deduced. 

VI. Discussion 

We have seen already in Theorem [T] how the regions of 
attraction constructed in Lemma [T] can lead to global conver¬ 
gence results. Here we discuss how our approaches connect 
with other results in the literature and provide also additional 
examples where delay-independent stability can be deduced. 

A- Connections with Lemmas \^and^ 

Lemmas |3] and H are important classical results that have 
been used in a number of interesting studies to prove delay- 
independent global asymptotic stability for linear, homoge¬ 
neous, and subhomogeneous monotone systems. These results 
are based on the fact that in such systems global asymptotic 
stability of the undelayed system implies the existence of arbi¬ 
trarily large points where the vector field is negative. Lemma [T] 
shows that such points can be determined under weaker system 
assumptions. Even though these points (denoted as Q can 
in general be difficult to characterize algorithmically, their 
significance lies in the fact that they allow to prove delay- 
independent global asymptotic stability for broader classes of 
systems. 
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The following corollary and its proof make the connection 
with Lemmas [3 and |4] more explicit. The corollary shows that 
when a point is known where the vector field is negative, 
then the delay-independent region of attraction obtained from 
Lemma [T] has a simple form, recovering arbitrarily closely 
that obtained from Lemmas [3] and |4l Note also that Lemma [T] 
additionally allows the delays to be time-varying. 

Corollary 3: If v is any point in R” such that f{v) + 
g^{vl) ^ 0 , then any solution of the delayed sys¬ 
tem (Ell with initial condition xtg € C([—r, 0], [0,u)) satisfies 
limt^oo x{t) = 0 . 

Proof: Define y{t) to be the trajectory of the undelayed 
system (|4]i through y(0) = v. Then the given vector field 
property implies that, for all sufficiently small > 0 , y(f^) 
y{0) holds and the function defined in (| 6 |l satisfies hi{tP) = 
for all i. Therefore, it follows from ® that the point ( can 
be brought arbitrarily close to v, and so [ 0 ,u) is a region of 
attraction obtained from Lemma [T] ■ 

Remark 6: Corollary [3 is also true if the inequality for the 
vector field is not strict, i.e. f{v)+g^{v'l) < 0 , and the point v 
can also be included in the estimate of the region of attraction. 
This follows directly from Step 3 of the proof of Lemma [T] 
by replacing C, with v. 


B. Subhomogeneous systems 

We now indicate how the regions of attraction constructed in 
Lemma [T] can also recover known delay-independent stability 
results for monotone systems that are also subhomogeneous. 
To see this, let us first recall the definition of subhomogeneity, 
as used in CD. 

Definition 1: The functions / and p* are said to be subho¬ 
mogeneous of degree a > 0 if they satisfy f{Xx) < A“/(x) 
and g*{X4>) < A“(;*(</)) for all A > 1, all < > 0, all x € R", 
and all G C([—r, 0], R"). 

Corollary 4: Suppose that / and p* are both continuously 
differentiable for all fixed t and subhomogeneous of any 
degree a > 0 and that the origin is known to be globally 
asymptotically stable for the undelayed system (|4]i- Then 
the origin is globally asymptotically stable for the delayed 
system ©. 

Proof: Consider the time-invariant bounding system (fT4l i. 
This is then positive, monotone, and subhomogeneous, and 
the origin of the corresponding undelayed system is known to 
be globally asymptotically stable, whence we recall from the 
proof of Theorem 3.2 in CD that given any xq > 0 , there ex¬ 
ists V ^ xq such that f{v)+g{vt) <C 0. Therefore, it follows 
by definition of g given in CD that f{v) + <C 0 , and 

so Corollary [3 guarantees that the region of attraction of the 
origin of the delayed system © contains the region [ 0 ,u). 
Allowing xo —oo, u can then be made arbitrarily large, 
meaning that the convergence region obtained from Lemma Cl 
is the entire positive orthant R". ■ 

The proofs of Corollaries |3] and 0] show that subhomoge¬ 
neous monotone systems admit undelayed trajectories through 
which Lemma Cl can be used to deduce delay-independent 
global asymptotically stability (with time-varying delays), 
recovering Corollary 3.3 from CD- It should be noted, though. 


as was shown in CD, that subhomogeneity is a sufficiently 
strong property that leads to delay-independent stability also 
for classes of systems that are not monotone. 


C. Global asymptotic stability via Gronwall’s Lemma 

We now give an example where the regions of attrac¬ 
tion quantified in Lemma Cl can be used to deduce delay- 
independent global asymptotic stability in a system where 
Theorem Cl does not appljo and which is also not subhomoge¬ 
neous. In particular we consider a delayed form of the second 
example in CD, with delays Ti{t),T 2 (t) G [ 0 ,r], 


dxi 

dt 


it) = -xiif) + 


X2{t - ri(f))^ 

X2{t - Ti(f))2 -p 1 

2 X 2 (f)^ 


^(t) =Xi(t-T2{t)) - , „ 

dt X2\ty -\-1 


( 18 ) 


The equations ( fTSl l define a monotone system for initial condi¬ 
tions taking values in R”, and the origin of the corre^onding 
undelayed system is globally asymptotically stablqj as dis¬ 
cussed in CD- For any k > 1, consider the initial condition 
y( 0 ) = in Lemma Cl It can be showr@ that 

the points can be chosen to satisfy C, S> (1 ,k — 1 )^. 
Therefore, by choosing k arbitrarily large. Lemma Cl gives 
delay-independent convergence on a region of attraction con¬ 
taining the set {x: Xi < 1}. Additionally, it can be deduced 
from an application of GronwalTs Lemma that all trajectories 
of the delayed system must reach this set in finite tim^B- 
Consequently, all solutions eventually enter a region on which 
Lemma Cl guarantees convergence. Therefore, the origin in 
the system CD* is globally asymptotically stable for arbitrary 
bounded time-varying delays. 


VII. Conclusions 

Within this paper we have considered the question of when 
it is possible to guarantee delay-independent stability for non¬ 
linear monotone systems of delay differential equations. We 

^Note that this example also does not satisfy [m Assumption A], which 
was used in (H to deduce global asymptotic stability for planar systems. 
Additionally, because both components in GD cannot simultaneously tend 
to oo in negative time. Theorem ^ also cannot be applied for this system. 

^Recall that the stability definitions considered are restricted to the non¬ 
negative orthant RJ. 

^Consider the undelayed system corresponding to Gl)., Whenever xi > 1, 


have < 0. Also, > 0 when xi € ( 


^+ 1 ’ 


i+i 


and ■ 


< 0 


when xi € (l, 




). Furthermore, for any 77 G (O, — l) and ^ > 




dx2 

dxi 




£ > 1, It holds that 

^ ’ \ ax-i\ I ^ I I ^ 

The above imply, for the trajectory y{t) through the given initial conditions, 
that the initial increase in y 2 until y\ = 2772/(^2 + 1) is smaller than its 
subsequent decrease until y\ = 1. Therefore, there exists fP > 0 for which 
(!,«: — 1)^ <C ‘C y{0), whence the point ^ defined in can be 

chosen such that ^ (1, ac — 1)^. 

^To see this, observe that whenever xi{t) > we have 
implying that it takes only finite time to reach the set {x: xi < |}. For any 
given trajectory of the delayed system, call this time r. Thereafter, the trajec¬ 
tory always remains within this set and hence (t) < | ~ which 

implies that X 2 {t) < b for all t > r, for some bounding constant b. Thus, 
we have X2{t) < B := max{sup^^[Q a72(f), < 00 . Consequently, 

< —xi{t) + -gfqrY all f > 0 and therefore, since < 1, 

Gronwall’s Lemma establishes the desired property. 
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showed that, when the undelayed system admits a convergent 
solution that is unbounded in all components in negative 
time, then this is sufficient to deduce that the corresponding 
delayed system will be globally asymptotically stable under 
arbitrary bounded delays, which can be heterogeneous and 
time-varying. This result followed from a more general result 
in which we showed that positive convergent trajectories of 
the undelayed system explicitly determine delay-independent 
regions on which all solutions of the corresponding delayed 
system must be convergent. These results demonstrate that it is 
possible to use the quasimonotonicity property alone to infer 
information about the behavior of a system under arbitrary 
bounded delays directly from knowledge of the asymptotic 
behavior of the corresponding undelayed system. 

Appendix A 

Quasimonotonicity and comparison properties 

The monograph provides an excellent reference on the 
general theory of monotone systems. Chapter 5 focuses on 
monotone systems of delay differential equations, analyzing 
their properties in detail and eventually proving a powerful 
result of generic convergence for autonomous systems under 
assumptions of boundedness, continuous differentiability, and 
irreducibility. Several of the preliminary results, however, are 
stated for general monotone delayed systems, making them 
useful here. We will now review two useful comparison results 
that are invoked within the proofs of our main results in 
Section lY] 

We shall consider the general delay differential equation 
^(f) =/(f,xt), (19) 

with / : A X C([—r,0],n) —)■ R" continuous on open 
subsets A C R and fl C R". We say that / satisfies the 
quasimonotonicity property if 

(/) < i/j and ^ ( 20 ) 

for all f € A and all & C([—r,0],n). Given Iq € A and 
(j) € C([—r, 0 ],n), we let f) denote the maximally 

defined solution that satisfies ( fTOl l for all t > tg and passes 
through Xtg = (fi¬ 
ll is well-known that systems that satisfy the property of 
quasimonotonicity are monotone, i.e. if two initial conditions 
(f>,ip £ C([—r, 0 ],r 2 ) satisfy <f> < if then x{t,to,(l), f) < 
x{t, to, ip, /) for t > to- This is stated in the following lemma 
(Theorem 5.1.1 in || 6 ]) which is the first result we recall. 
Note that this lemma also allows trajectories of systems with 
different vector fields to be compared. 

Lemma 3: Let f,g:Ax C{[—r, 0], fl) —>• R" be continuous 
in t and locally Lipschitz in cp, and assume that either f or g 
satisfies (Eil. Assume also that f{t, cp) < g{t, (p) for all f C A 
and all cp £ C([—r, 0], fl). If fo C A and (p,ip £ C([—r, 0], fl) 
satisfy (p < ip, then x{t,to,(p, f) < x{t,to,ip,g) holds for all 
t > to for which both are defined. 

The second result guarantees monotonic convergence of any 
bounded trajectory of an autonomous system at whose initial 
condition the vector field is either nonnegative or nonpositive 
in all components. This is Corollary 5.2.2 in [| 6 ]. 


Lemma 4: Let / : C([—r, 0],n) —?• R" be time-invariant, 
locally Lipschitz, and satisfy (l20l l. If w C 17 is such that 
f{vl) > 0 (f{vt) < 0), then x{t,to,vl, f) is nondecreasing 
(nonincreasing) in f > 7o- If the positive orbit of hi has 
compact closure in 17, then there exists k > v (k < v) such 
that x{t, to,vl, /) —>■ fc as 7 ^ cx). 

Appendix B 
Nonpositive systems 

The approach within the main paper focused, for simplicity 
of presentation, on monotone systems with an equilibrium at 
the origin. We now demonstrate that the analysis in Section 
|IV]can also be applied to yield delay-independent convergence 
guarantees in general monotone systems of the form of (|4|i 
and (O admitting an arbitrary equilibrium point x* £ R". We 
thus use within this section the standard definitions for stability 
as stated, for example, in Definition 5.1.1 in ifTSl . 

The monotonicity properties O and © imply that the 
orthants x* + R" and x* — R" must be positively invariant for 
both systems dUl and ©. It is therefore natural that coordinate 
transformations can enable Theorem [T] to be extended to the 
following more general result. 

Corollary 5: Suppose that the undelayed system (Hll ad¬ 
mits solutions y{t) and y{t) satisfying lim(_>._oo ^^(7) = 
oo, limt^ooyj(7) = X*, limt^_oo y.(7) = -oo, and 

lim(_>.oo 2/j(f) = X* for components i. Then x* is globally 
asymptotically stable for the delayed system ©. 

Analogously to Section HYl this represents a limiting case 
of a more general technical result which extends Lemma [U to 
quantify delay-independent regions of attraction. Note that the 
points (p, ( mentioned in the corollary are defined within the 
proof of the corollary in a way analogous to the point ( in 
Lemma [T] 

Corollary 6: Suppose that the undelayed system (HI admits 
solutions y{t) and y(7) satisfying ^^(0) > x*, limt_>oo (f) = 
X*, y-{0) < X* , and limt_).oo 2/^(7) = x* for all components 
i. Then there exist points C, > x* and C, < x*, defined in 
terms of the trajectories y{t) and y(t) respectively, such that 
any solution of the delayed system © with initial condition 
xto G C([-r, 0], [C,C]) satisfies limt_>oo x(7) = x*. 

If the vector field is additionally continuously differentiable. 
Corollary | 6 ] again extends as follows to enable stronger con¬ 
clusions to be drawn. 

Corollary 7: Suppose that / and g* are both continuously 
differentiable for all fixed 7 and that the undelayed sys¬ 
tem @ admits solutions y{t) and y{t) satisfying y^fO) > x*, 
limt_>oo = X*, y^{0) < x*, and limt^^y^ft) = x* for 
all components i. Then the points (f, ( in Corollary | 6 ] can be 
chosen such that C 2 > x* and <C x*, and the point x* is 
asymptotically stable for the delayed system ©. 

Corollary |7] then also immediately implies the following 
simple corollary. 

Corollary 8: Suppose that / and g* are both continuously 
differentiable for all fixed 7 and that the point x* is asymp¬ 
totically stable for the undelayed system I©. Then x* is also 
asymptotically stable for the delayed system ©. 


We now show how the results stated above can be derived 
from the theory established in Sections |IV] and IVl 

Proof of Corollary^: Make the coordinate change u = 
X — X* 'm ^ and Q. This gives 

+ u{t)) + g\x*t + u(t)t) 

= +9\u{f)t) ( 21 ) 

and 

= f{x* + u{t)) + g\x*t + ut) 
at 

= + 9\ut). (22) 

The systems (I2TI1 and (l2^ then have an equilibrium at the 
origin and, moreover, their nonlinearities / and g* clearly 
satisfy the monotonicity properties ([T]i and (O. Finally, using 
the given trajectory y{t) of (| 4 |, we see that (|2TI) admits a 
solution v{t) defined by v(f) = y{t) — x* for all f > 0 
satisfying t 7 i( 0 ) > 0 and limj^oo Ui(f) = 0 for all components 
i. The analysis of Lemma [T] can thus be applied, yielding a 
point 7 > 0 dehned as in ([8]l such that the origin of (|2^ 
is asymptotically stable with region of attraction contain¬ 
ing [0,7]. 

Analogous analysis with the positive monotone systems 

flU 

-{t) = -f{x*-u{t))-g\xn-u{t)l) 

= Ruf)) + ( 23 ) 

and 

nil 

fjfi.t) = -f(x* - u{t)) - g\x*l - Ut) 

= fiRt))+Riut), ( 24 ) 

obtained under the coordinate change u = x* — x, gives a 
point 7 > 0 such that (l 24 ll has asymptotically stable origin on 
the region of attraction [0,7]. 

Combining these analyses and denoting ( = x* +j and 
( = X* — y shows that the region of attraction of the equilib¬ 
rium point X* of the original delayed system Q contains the 
union [C, a:*] U Finally, let z(f) and z(f) respectively 

denote the trajectories of Q with initial conditions (^1 and (^1 
and let also C ^ a;to < C- Then the fact that z(f) and z{t) 
converge to x* in conjunction with Lemma 0 guarantees that 
limt_).oo x{t) = X* also. Therefore, the region of attraction 
must in fact contain [C,C]- ® 

Proof of Corollary [ 7 } Analogously to the proof of 
Corollary [H we invoke ||6l Remark 3 . 1 . 2 ] to deduce that the 
points in Corollary |6] satisfy C ^ a:* ^ C- ■ 

Proof of Corollary]^ As in the proof of Corollary | 2 ] the 
asymptotic stability of (01 implies the existence of e > 0 such 


that limj^oo y{t) = 0 for solution y{t) of (01 with initial 
condition ?/(0) G {x*,x* -he). Analogously, there exists e > 0 
such that limt^oo2/(f) = 0 whenever y{0) G {x* — e,x*). 
Corollary [T] can then be applied to deduce the result. ■ 
Proof of Corollary |5].' Analogously to the proof of 
Corollary |6l the trajectories v{t) = y{f) — x* and yit) = 
x*—y(t) satisfy the asymptotic properties required in Theorem 
[T]for the positive monotone undelayed systems (1211 1 and ( l23l l. 
respectively. Therefore, applying Theorem [T] and invoking 
Lemma |3] as before immediately implies the desired result. 


References 

[1] A. Kolmogorov, “Sulla teoria di VolteiTa della lotta per I’Esistenza,” 
Giorno Istituto Italiano degli Attuari, vol. 7, pp. 74-80, 1936. 

12] C. Altafini, “Consensus problems on net\vorks with antagonistic inter¬ 
actions,” IEEE Trans. Autom. Control, vol. 58, no. 4, pp. 935-946, 
Apr. 2013. 

13] R. D. Yates, “A framework for uplink power control in cellulai' radio 
systems,” IEEE J. Sel. Areas Commun., vol. 49, no. 7, pp. 1341-1347, 
Sep. 1995. 

14] M. W. Hirsch, “Systems of differential equations which are competitive 
or cooperative,” SIAM Journal of Mathematical Analysis, vol. 13, no. 2, 
pp. 167-179, Feb. 1982. 

15] M. W. Hirsch and H. Smith, “Monotone dynamical systems,” Handbook 
of Differential Equations: Ordinary Differential Equations, vol. 2, pp. 
239-357, Sep. 2005. 

16] H. L. Smith, Monotone Dynamical Systems: an Introduction to the The¬ 
ory of Competitive and Cooperative Systems. American Mathematical 
Society, 1995. 

17] L. Farina and S. Rinaldi. Positive Linear Systems, Theory and Applica¬ 
tions. John Wiley & Sons, 2000. 

18] X. Liu, W. Yu, and L. Wang, “Stability analysis for continuous-time 
positive systems with time-varying delays,” IEEE Trans. Autom. Control, 
vol. 55, no. 4, pp. 1024-1028, Apr. 2010. 

19] O. Mason and M. Verwoerd, “Observations on the stability properties 
of cooperative systems,” Systems and Control Letters, vol. 58, no. 6, pp. 
461^67, Jun. 2009. 

[10] V. S. Bokharaie, O. Mason, and M. Verwoerd, “D-stability and delay- 
independent stability of homogeneous cooperative systems,” IEEE Trans. 
Autom. Control, vol. 55, no. 12, pp. 2882-2885, Dec. 2010. 

[11] V. S. Bokharaie and O. Mason, “On delay-independent stability of a class 
of nonlinear positive time-delay systems,” IEEE Trans. Autom. Control, 
vol. 59, no. 7, pp. 1974-1977, Jul. 2014. 

[12] V. S. Bokharaie, O. Mason, and M. VeiAvoerd, “Correction to D-stability 
and delay-independent stability of homogeneous cooperative systems,” 
IEEE Trans. Autom. Control, vol. 56, no. 6, pp. 1489, Jun. 2011. 

[13] J. K. Hale and S. M. Lunel, Introduction to Functional Differential 
Equations. Springer-Verlag New York, 1993. 

[14] A. Rantzer, B. Riiffer, and G. Dirr, “Separable Lyapunov functions for 
monotone systems,” Proc. 52"^ IEEE Conf. Decision Control, pp. 4590- 
4594, Dec. 2013. 

[15] J. D. Murray, Mathematical Biology I: An Introduction. Springer, 2002. 

[16] M. Haque, S. Sarwardi, S. Preston, and E. Venturino, “Effect of delay in 
a Lotka-Volterra type predator-prey model with a transmissible disease 
in the predator species,” Mathematical Biosciences, vol. 234, no. 1, pp. 
47-54, Nov. 2011. 

[17] J. Hofbauer and J. W. H. So, “Multiple limit cycles for three dimensional 
Lotka-Volterra equations,” Applied Mathematics Letters, vol. 7, no. 6, 
pp. 65-70, 1994. 

[18] E. Devane and 1. Lestas, “Delay-independent asymptotic stability in 
monotone systems,” Proc. Amer. Control Conf, Jun. 2015. 


